Simulation of Coherent Synchrotron Radiation Emission from Rotating Relativistic 

Electron Layers 
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The electromagnetic radiation of rotating relativistic electron layers is studied numerically using 
particle-in-cell simulation. This system is general enough to be applied to problems involving the 
coherent synchrotron radiation (CSR) instability in bunch compressors and possibly in radio pulsars 
where the CSR instability may be used to explain the high brightness temperature and the observed 
spectra. The results of the simulation confirm all relevant scaling properties predicted by theoretical 
models and suggest that their range of validity is bigger than conservative estimates indicate. This 
is important for radio pulsars where the allowed parameter range is unknown. 



Extremely high brightness temperatures encountered 
for example in the radio emission of pulsars cannot be ex- 
plained by an incoherent radiation mechanism for which 
the radiated power is proportional to the number of ra- 
diating charges. However, if all outgoing waves interfere 
constructively the total power scales as the square of that 
number. In future particle accelerators the wavelength 
of the emitted radiation can easily reach the length of 
the accelerated bunches such that bunches can modulate 
themselves by means of their own radiation. The possi- 
ble energy loss caused by CSR would be undesirable for 
the operation of such accelerators. 

The aim of the present work is to test the results ob- 
tained from the model developed by Schmekel, Lovelace 
and Wasserman Q numerically and probe its range of 
validity. The model predicts growth rates, saturation 
amplitudes and the emitted power for an initial pertur- 
bation of a rotating relativistic cylinder of charged parti- 
cles in an external magnetic field. The linear instability 
is closely related to the instability found by Goldreich 
and Keeley Q. More recently CSR has been investi- 
ated in the particle accelerator physics community, e.g. 



In [l| the distribution function is chosen such that 
the azimuthal canonical angular momentum is fixed and 
the energy drops off exponentially. For the time be- 
ing all fields and distribution functions are assumed to 
be uniform in the z-direction. The equilibrium can be 
parametrized by the radius r of the cylinder, the relative 



energy spread 5E/E 



of the particles, the Lorentz 



factor 7 and the dimcnsionlcss field reversal parameter 



(i) 



where no is the central number density, the az- 
imuthal velocity, B^ xt is the external magnetic field 
in the z-direction in cgs units and c — 1. Consid- 
ering all perturbed quantities to have the dependence 
exp (imcj) — iut + k r r) one obtains for m 2 / 3 v t h <C 1, 
krVth.ro <C 1 and \Au>\ <C 7 -2 the dispersion relation 



and the growth rate 



l/2 m 2/3, 



1.083C i/2 m 
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where Aw = (u> — m<f>) / '(mcf)) , H m = J m + iY m with the 
Bessel functions J m and Y m , (f> is the angular velocity, 
the dot denotes the time derivative, and 



Fo = — 
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The radiated power is given by 



1 



-mC,J m {ujro)H m (ujro)~i 2 (Au)) 2 F 



(2) 



P m «3.71xl0 1 Vm- 3 -f^ > ) ( 5 ) 
ro V <t> J s 

Here the length of the cylindrical layer is denoted by 
L. The quoted expression was derived by approximating 
Bessel functions with Airy functions. It was shown that 
this approximation is valid under the conditions m 3> 1 
and m 3> £ 3 / 2 7 3 = m\. The used dispersion relation 
was derived under the assumption |AcD| <§; -f~ 2 (which is 
approximately equivalent tom>mi), but this condition 
may be too strict. There is evidence that this is indeed 
the case and one objective of this article is to confirm this. 
Especially for the study of radio pulsars bigger growth 
rates may be necessary to explain the observed brightness 
temperature. 

For the simulation the software package OOPIC was 
used. OOPIC is a relativistic two-dimensional particle- 
in-cell code which supports both plain (x, y)-geometries 
and cylindrical (r, z)-geometries. Since the interesting 
dynamics takes places in the azimuthal direction one can 
only simulate a thin ring (instead of a cylinder) in the 
(x, ?/)-mode. Loading the initial circular particle distri- 
bution in the (x, y)-mode required modifying the source 
code (files load.cpp, diagn.cpp and c_utils.c) to allow the 
program to handle circular particle distributions. Some 
minor modifications were necessary in order to compile 
XOOPIC-2.5.1 with gec 3.2.2 and the compiler compiler 
bison 1.28 under SunOS 5.9. The built-in function parser 
was extended to support elliptic integrals. 
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Since a thin ring of particles is simulated instead of 
a cylinder thereof all fields and charges were divided by 
the length L of the cylinder whereas the electron mass 
needs to be divided by L 2 . L = 10r is chosen unless 
noted otherwise. The electric and magnetic self-fields 
for a thin ring equilibrium differ from what was used in 
the model. The fields can be found in 8]. It is ensured 
that OOPIC uses these self-fields before the perturbation 
starts to build up. As it turns out choosing the correct 
self-fields is not too crucial. Leaving them out the system 
will build them up itself. Once the self-fields are created 
the system shows no difference in behavior. The absence 
of the self-fields in the dispersion relation might help to 
understand this feature. As in 1] a Gaussian number 
density profile with RMS width vth was chosen for the 
initial distribution. 5000 macro particles were tracked on 
a grid with resolution 512 x 512 unless noted otherwise. 
Once an energy for a particle has been chosen it is placed 
at the equilibrium radius rg = m7c(eB) _1 , i.e. neglect- 
ing betatron oscillations particles on the same orbit have 
the same energy. This fixes the azimuthal component of 
the canonical angular momentum. The system can pick 
up transverse motion quickly. The grid represents a rect- 
angular region 40m x 40m big where the ring with radius 
ro = 10m is centered. 

In Fig. ^ the initial particle distribution (gray) and 
the particle distribution after 23ns are shown. The pa- 
rameters are ( = 0.010, 7 = 30, and v t h — 0.002. Qual- 
itatively, a bunching of the particle distribution can be 
observed. An enlargement of a small section of Fig. Q]is 
also shown in the same figure. In Fig. [21 the bunching is 
shown for successively higher energy spreads. With in- 
creasing energy spread the bunches become fuzzier and 
the clean gaps between bunches that can be observed 
for small energy spreads are populated with "stray par- 
ticles". This suggests that it may be harder to achieve 
complete coherence for larger values of v t h- The deco- 
herence due to the non-zero width of the particle beam 
is investigated quantitatively later in the paper. These 
qualitative features are independent of 7. 

Also note that during the evolution of the circular 
charge distribution both the radius and the width of the 
ring increase slightly. The former is due to a.) particles 
losing energy and b.) the perturbed magnetic field chang- 
ing significantly. It tends to decrease for small energy 
spreads and increase for larger energy spreads. Starting 
with a larger radius the radius increases even further, i.e. 
this is not a relaxation from a "false" to a true equilib- 
rium. Since the non-zero mesh size imposes an upper 
limit on the azimuthal mode number m which can be- 
come unstable, it is expected that the distance between 
bunches decreases as the resolution increases. This is in- 
deed the case. For larger energy spreads the bunching of 
the distribution becomes hardly visible, but it still can be 
observed in the z-component of the magnetic field (Fig. 

EJ). 
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FIG. 1: Initial particle distribution (gray) and the same dis- 
tribution after 23ns have elapsed. Parameters: £ = 0.010, 
7 = 30 and v th = 0.002 



•ST 

".if: 



„• - » •. V 



FIG. 2: Particle distribution (7 = 30, ( = 0.01) after 23ns for 
Vth = 0.002, v th = 0.008, v th = 0.015 and v th = 0.033 (from 
left to right). The four figures depict an approximately 6.5m 
long segment. 



The bunches are slightly tilted and may be connected 
by a very thin inner ring of particles for sufficiently high 
beam currents. For these reasons it is not possible to 
Fourier transform the charge perturbations in order to 
compute the growth rates for each value of m. Since the 
resolutions used were low the range of m values is re- 
stricted. Therefore, only the radiated power is computed 
which can be obtained easily. Estimates of the growth 
rate are two orders of magnitude higher than what would 
be expected. A possible explanation is that the ratio be- 
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FIG. 4: Loss of kinetic energy in W vs. 7 for £ = 0.01 and 
v th = 0.002. 
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FIG. 3: z-component of the magnetic field (self-field plus 
perturbation without external magnetic field) after 23ns for 
C = 0.010, 7 = 30 and v th = 0.025. The size of the area 
depicted is 40m x 40m. 



tween the saturation amplitude and the electric self field 



E. 



self 



1 /3(w)(m) 



(0) 



is typically in the order of 10~ 3 for the given sample cases 
which is rather small. The initial perturbations due to 
discreteness, numerical noise etc. are usually in the same 
order of magnitude. Therefore, one cannot expect to 
see the regime covered by the linearized Vlasov equation. 
This is another reason for focusing entirely on the emitted 
power. 

In Fig. 0|and Fig. [3]the radiated power determined by 
measuring the kinetic energy loss of the electron cloud 
after approximately 2.36 ns is plotted as a function of 7 
and C, respectively. After 0.24 ns the perturbations have 
saturated and the emitted power is fairly constant. A 
quadratic dependence can be established, i.e. the first 
two relevant scalings expected from the analytical model 
are recovered. 

The simulation has been repeated at lower (256 x 256) 
and higher (1024 x 1024) resolution. No significant effect 
could be observed. This is consistent with the model 
which predicts that most power is emitted by modes with 
low m. Also, decreasing the stepsize dt to 2.5 ps and 
increasing the number of macro particles to 50000 has a 
negligible effect. 
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FIG. 5: Loss of kinetic energy in W vs. £ for 7 
vth = 0.025. The solid line shows the best fit. 
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Finally, the effect of the energy spread is investigated. 
With increasing vth the power decreases which is due to 
the decoherence described by the factor Fq defined in Eq. 
0] The results are plotted in Fig. HJ for the parameters 
C = 0.01, 7 = 30 and C = 0.005, 7 = 10, respectively, 
and L = ro- In the former case mi is 27 and in the latter 
case it is 0.4. Eq. becomes 



P 



,14 erg ijv^, 
s r Q ^ 



3.71 x io«:=° _y > m[mJ m {Lor )H m {u;r )Y (1) 
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FIG. 6: Total power radiated as obtained from OOPIC for the 
parameters C = 0.01, 7 = 30 (solid) and £ = 0.005, 7 = 10 
(dashed). The dash-dotted line is proportional to X^m=i 

Despite m\ being much larger than 1 for the first set 
of parameters the slopes in Fig. H3 match exactly only 
if the summation starts at m — 1. This suggests that 
modes with m < mi do radiate and can be described 
by the same dispersion relation. Since the power scales 
as m~ 5 / 3 these modes may actually be very important 
for computing the total energy loss. Note that while the 
simulation suggests P oc L 2 Eq. [S] (which was derived 
under the assumption L > ro) gives P oc L. A 2D sim- 
ulation cannot explain how the radiation from different 
axial positions on the cylinder interacts. In the thin ring 
case doubling L doubles the number of particles N and 
therefore quadruples P. Fortunately, as can be seen in 
the derivation of Eq. |3]in Q the C, Vth and 7 dependent 
part of P is independent of L. In Fig. EJthe overall factor 
matches if ro = 100-L whereas the growth rate for a per- 
turbation of a cylinder and a thin ring coincide for tq = L 
0. Also note that Fig. Q]suggests k r vthfo ~ 1, whereas 
Eq. 13 was derived under the assumption k r v t hTQ <C 1. 



The particle in cell code OOPIC was used to simulate 
the evolution of density perturbations in a thin ring of 
charged particles which move in relativistic almost cir- 
cular motion in an external magnetic field. The results 
were compared with the model in ij . Comparisons of the 
simulation with the model shows approximate agreement 
with the main predicted scaling relations. In particular 
the bunching effect could be observed very clearly and 
the emitted power is proportional to the square on the 
number density which implies coherent radiation. The 
dependence on the energy spread can be recovered ex- 
actly assuming all modes contribute to the observed en- 
ergy loss suggesting that the model may apply even if 
m < mi. 
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